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Abstract 

The dynamics of a coupled system comprised of a two-level atom and cavity field assisted by 
continuous external classical field (driving Jaynes-Cummings model) is studied. When the initial 
field is prepared in a coherent state, the dynamics strongly depends on the algebraic sum of both 
fields. If this a sum is zero (the compensative case) in the system only the vacuum Rabi oscillations 
occur. The results with the dissipation and external field detuning from the cavity field are also 
discussed. 
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I. INTRODUCTION 



The ability to create, manipulate, and characterize quantum states is becoming an in- 
creasingly important area of physics research, with implications for such areas of technology 
as quantum computing, quantum cryptography, and communications, see jlj], Q], Q], 
Most of research in quantum nonlocality and quantum information is based on entangle- 
ment of two-level particles. One of the most interesting aspects of its dynamics is the 
entanglement between atom and field states. This essentially quantum mechanical property 
with no classical analog is characterized by the impossibility of completely specifying the 
state of the global system through the complete knowledge of the individual subsystem's 
dynamics. 

The Jaynes-Cummings model[5] (JCM) for the interaction between a two-state atom and 
a single mode of the electromagnetic field holds a central place in description of such inter- 
action and provides important insight into the dynamical behavior of atom and quantized 
field. In driving JCM the cavity field and driving field start to interact, which provide a 
possibility to study directly the field dynamics at joint interaction with a two-level atom. 
Recently, it was shown that the effective coupling between an atom and a single cavity field 
mode in JCM '.( driving JCM) can be drastically modified in the presence of a strong external 
driving field fj,[zj], 0. The important line of this direction is to use microcavities and 
microspheres for changing the features of atom-field interaction as a result of placing atom 
or quantum dots into a microcavity (seeflfll]. jll]. l^). 

The driven Jaynes-Cummings model for cases where the cavity and external driving 
field are close to or in resonance with the atom, has been studied by several authors. In 
Ref. [3] studied the Stark splittings in the quasienergies of the dressed states resulting from 
the presence of the driving field in the case where both fields are resonant with the atom. 
Authors studied the effect of the external field on the Rabi oscillations in the case where 
the cavity field is resonant with the atom and where the external field is both resonant 
and nonresonant. In Ref. [lfij] studied a similar model where the external field was taken to 
be quan tized. Much attention was given to the limit of high-intensity of driving field. In 



Ref. 



16 1 have studied the JCM with an external resonant driving field and have shown that 



the collapses and revivals of the mean photon number occur over a much longer time scale 

n 

than the revival time of the Rabi oscillations for the atomic inversion. Author 17] studied 
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the interaction of an atom with both a quantized cavity field and an external classical 
driving field, in the regime where the atom and fields are highly detuned. He has shown how 
dispersive interaction can be used to generate coherent states of the cavity field and various 
forms of superpositions of macroscopically distinct states. 

The main goal of the present work is the calculation of the dynamics of coupled atom and 
field assisted by continuous external pumping in case when the initial field is prepared in a 
coherent state. The main result is the following: Starting with a field's mode in a coherent 
state and with the atom in its upper state, the dynamics strongly depends on the algebraic 
sum of amplitudes of initial cavity field and the external field. If such sum is close to a zero 
(the compensative case), in system only the vacuum Rabi oscillations occur. 

This Letter is organized as follows. In Section 2 we discuss the motion equations for two- 
level atom coupled to the field in cavity with the assistance of continue pumping classical 
field. Section 3 presents the results of numerical study for the dynamics of the atom and 
field subsystems for dissipative case by the technique of the master equation. The behavior 
of entropy and Fourier spectrum of oscillations is studied also. In last Section, we discuss 
and summarize our results. 

II. BASIC EQUATIONS 

Consider a two-state atom, driven by a classical external field (l/2)E e exp(iu e i) + c.c, 
and coupled to a cavity mode of the quantized electromagnetic field. The Hamiltonian for 
the atom-cavity system (assuming h — 1) in rotating- wave approximation (RWA) is given 
by 

H = -u a 3 + u a a + a + g [a~a + + a + a] + - [Ea+e^ 1 + ^Ve - ***] , (1) 

where ojq is the transition atom frequency, u a is the cavity frequency, g is the coupling 
constant between the atom and the cavity field mode, £ is proportional to the coupling 
constant between the atom and the external classical field of frequency u e and the amplitude 
of that field, a + and a are the creation and annihilation operators for the cavity mode 
[a, a + ] = 1. In general uj q , u a , and uj e are different. To remove the time dependence in H, 
we use the operator exp[— icu e t(a 3 + a + a)] to transform to a frame rotating at the frequency 
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u e . The Hamiltonian in the rotating frame (the interaction picture) is then 



Hi = -(T3+9 [o-a + e~ m + a + ae i5t ] + ~ [£a+ + £*a~] 



(2) 



2 J 2 

where A = — u e , 5 = u) a — u e , a ± = (cr x ± ia v ) /2, cr z , cr y , cr 2 are Pauli matrices. First 
we consider the resonant case when A = and 5 = 0. Case of non-zero detuning 5 ^ is 
discussed in the second part. The resonant Hamiltonian (ou a = ujq = u e ) in the interaction 
picture has the form 



Hi, 



g [<j-a + + a + a] + - [£a + + £*a~] . (3) 

Further we use the following dimensionless variables r = u Q t, g/^o — » g, £/^o — ► £■ If £ = 
the Eq.fjnj) describes the standard JCM, case £ ^ corresponds to driving JCM. To obtain 
the solution to Eq.Q we introduce the displacement operator -D( 7 ) = exp{7<2 + — 7*a}, 
7 = £/2g, which allows us to rewrite Eq.Q in the form 



H ir = gD + (j)((T + a + a+O^W, 



(4) 



where identity D + ('~f)(a + , a)£)( 7 ) = (a + +7*, a + 7) is used. Establishing in (J3J) the Hamil- 
tonian D{ r ))H ir D + {^j) = g(a + a + a + a~) and the state vector ip) = -D( 7 ) we obtain the 



well-known the Schrodinger equation of the standard JCM in form 







Or 



[a a + a a~ 



(5) 



Now consider the case when the initial state of the field in cavity is a coherent state \a), 
with a = n 1 / 2 e~ m (n is the average number of photons in the field). Also we assume the 
atom is prepared in the excited state \e) (\g) is the ground state). The initial state vector 
\ip) = \e) \a) = \e) D(a) |0) allows us to write for Eq.(jSJ) the corresponding initial state 



vector 



ip ) as 



^ = D( 7 ) = \e) D(j)D(a) |0) = |e> | 7 ) , (6) 
where 7 = 7 + and overall factor exp(z Im( 7 a*)) is dropped. With Eq.(jnj) the solution to 
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the standard JCM Eq.(j3J) is given by 



^ (0) = ^C n (7){cos(eV^+T) |e) \n) - i sin^v^TT) \g) \n + 1)}, £ = <?r, (7) 



n=0 



where C n = C n {^) = exp(— |7| 2 /2)7 n /\/n[ are expansion coefficients for (7) state in the 
number representation \n). The Eq. (JZJ) allows us to write the solution to resonant driving 



ip ) in the following form 



JCM © |V(0> = D{~l) 



IV>(£)> = ^C„(7){cos(ev^Tl) |e> |- 7 ;n) -isin(eV^TT) |<7> |- 7 ;n+l», (8) 

n=0 

where |— j;n) = D(— 7) |n) is the displaced number state. Note the following. Formally 
quantities 8 and a have different physical meaning. Quantity S is used as a parameter of the 
driving field in the Hamiltonian while a is a factor of initial condition for field's mode 
in the cavity. However dependence of probability coefficients C n (^y) on the 7 = £/2g + a in 
(JHJ) shows a deep similarity of these quantities for the resonant case 6 = 0. The coherent 
field state has minimum uncertainty, and resembles the classical field as closely as quantum 
mechanics permits Q]. From the Eq.lfS]) the density operator p can be written as follows 

P = {HO) MO I = \e) (e\ U ee + \e) (g\ U eg + \g) (e\ U ge + \g) (g\ U gg , (9) 
where matrix elements C/y are given by 

00 

U ee = Y, C n C rnCOS^V^TT)cos^V^+l)\-T,rn)(-r,n\, (10) 

n, 771=0 
00 

U eg = —i C* n C m cos(£Vm + 1) sin(£vVi + 1) I -7; w) (-7; + 1| , 

71,771=0 
OO 

U eg = C* n C rn cos(£\/m + 1) sin(£\/n + 1) (-7; m + 1) (-7; n| , 

71,777=0 
OO 

= C*C m sin(£>/ m+T) sin(£\/ra + 1) |- 7 ;m + 1) (-7; n + 1| , 

71,771=0 

and are still operators with respect to the field. To describe the evolution of the atom (field) 
alone it is convenient to introduce the reduced density matrix 



P a(f) = Tr f{a) {p}, (11) 
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where the trace is over the field (atom) states. We have used the subscript a, f to denote the 
atom (field). Unlike the state vector, the density operator does not describe an individual 
system, but rather an ensemble of identically prepared atoms, see e.g. Ref.|3|. A condition 
for the ensemble to be in a pure state is that Tr{(p a, -Q 2 } = 1. In this case a state-vector 
description of each individual system of the ensemble is possible. On the other hand, for a 
two-level system, a maximally mixed ensemble corresponds to Tr{(p a '^) 2 } = 1/2. Due to 
the identity 

Tr f {\- T , m) (-7; n\} = Tr a {D + ( 1 ) \m) (n\ L>( 7 )} = Tr a {\m) (n\ L>( 7 )L> + ( 7 )} = 5 mn 
one can write p a in the following form 



P a = {J2 l^-i I' ( cos2 (^) \e) (e\ + sin 2 (£v^) \g) (g\)}+ (12) 

n=l 

00 

+ ff{^2\C n -i\ 2 fcos(£\/n + 1) sin(^v/n) |e) (g\ - sm(£Vn + 1) cos^y/n) \g) (e\)}. 

n=l 

From Eq.([12jl one can see that off-diagonal matrix elements p a eg = (e\ p a \g) and p a ge = (Pe 9 ) + 
are of the first order with respect to 7 = 7 + a, 7 = S/2g and therefore contain the 
information on the relative field's phase even in a weak field limit. The mean photon 
number (n) is given by 



(n) = (a + a)^, = Tr{a + ap^}. (13) 
Taking into account the identity 

Tr{a + a |— 7; m) (—7; n| } = Tr{D + Da + D + DaD + ('~f) \m) (n\ D(j)} = (n\ (a + — 7*)(a— 7) \m) 
and after minor algebra one may write ([TBI ) in the following form 

(n) = (a + a) f = A - B - B*, (14) 

where 



00 



A = |7| 2 + | 7 | 2 + ]T IC^f sm 2 (grV^), (15) 



n=l 
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1 _ 00 1 

~~ ra=l 



= Vri + 1 ± \fri, C n = exp(- |7| /2)—=, 7 = 7 + 0; = — + a. 



'n! 2p 

One can see, the quantum Rabi oscillations in appear as sum of the sinusoidal terms 
at incommensurable frequencies g^/n, gQ^i weighed by the probabilities |C n (7)| 2 . For the 
vacuum case a = equations (|14 p -([16 p are in agreement with Ref.[16]. One can see from 
equations (|15|). ()16|) the following. Since \C n \ 2 has maximum at n ~ |7| 2 (n > 1) one may 
emphasize a desired frequency in the spectrum the excited-state probability P + {t) = p a ee 
|T2|) (the probability that the atom is in the excited state) and (n) (fTl|) by a corresponding 
choice of complex quantity S. Simple calculation yields the number of corresponding Rabi 
frequency equals to E^ 2 — 1], E[x] is integer part of x. In other interesting case, the external 
field may be chosen as 

8 = -2ga. (17) 
In this case in equations (|T2|) . fT4j) 7 = 0, i? = 0, C n = 5 n \ , so P + and (n) reduce to the 



following simple form 



P = cos (gr) and (n) 



5 2 



+ sin 2 (sr). (18) 



2g 

In this case the two-level atom, which is initially in the excited state, undergoes the one- 
photon oscillations (radiating and absorption of a photon) and only the vacuum Rabi fre- 
quency peak is present in the frequency spectrum. Note the Rabi oscillations for standard 
JCM were directly observed in Ref. [2(|. 

The next experiment can be proposed on the basis of equations (|17j) -(|18| ) . One can vary 
both amplitude and phase of the external field 8 until only single vacuum Rabi frequency 
spectrum is observed. In this case condition a = —8/2g must hold, which provides the 
opportunity to measure parameters of the coherent state \a). In a compensative case (fTTj) 
these fields have equal amplitudes and frequency, but are shifted in phase by n. One can 
interpret this as follows: in the resonant case (5 = 0) the total field in the cavity may 
be written as a superposition of both fields in form (|o) ± |— a)) /V2. It generates the 
Schrodinger cat states since such field is periodically entangled with a two-level atom while 



7 



vacuum Rabi oscillations occur. Method of generating such states have been explored by a 
number of authors(see 21 1, 2^|, 3] and references therein). 

The above theory is valid in the case of very small dissipation (which describes the 
interaction of the atom and field subsystems with the environment) and zero detuning 5 = 0. 
But in the experiments, the damping of the cavity mode and the rate of spontaneous emission 
of the atom are not negligibly small. Thus, a cavity damping must be included in a treatment 
of driving JCM to compare it with experiments. An equation for density operator p(r) is 
required (master equation) because the loss of the coherence due to the reservoir changes 
any system pure states to mixed states. With cavity damping in effect we must solve master 
equation for the joint atom-field density operator p of the two-level atom coupled to the 
electromagnetic field in cavity. Such equation is given by 

^■ = -i[H i ,p]+C 1 p + C 2 p, (19) 

(XT 

where Hi is written in For simplicity we only consider the case 5^0 and A = 

0. Master equation ()19|) is more difficult to solve and numerical methods usually must 
be used. At interesting frequencies range the dissipation is written in Eq.(fTU|) both as 
mirror losses in the cavity that defines the mode of the electromagnetic field Cip, and 
as spontaneous emission from the atom £ 2 p. At Born-Markov approximation and zero 
temperature, these parts are written in equation (fT^j) as the following. One term is as follows 
C\p = 71 (2apa + — a + ap — pa + a), where 71 is the rate of single-photon losses. Another term 
£-2P = (72/2) (2apa + — a + ap — pcr + a) takes into account the spontaneous emission from the 
atom out of the sides of the cavity |23|]. In this case the atom is damped by spontaneous 
emission with damping rate 72 to modes other than the privileged cavity mode with frequency 
u a . We have solved the equation (fTHjl numerically using a truncated number states \n) and 
atom states |e), \g) basis. The algorithms for integration of such a system numerically can 



be found, e.g. in Ref . [24j . In general < n < 00, but for numerical calculations we have 
used < n < M . A finite base of the number states M was kept large enough so that the 
highest energy state is never populated significantly. Since we assume above that the field 
is initially in the coherent state |a), and the atom is in the upper state |e), then initially 
p = \a) (a I <S> |e) (e| must hold. Using numerically obtained joint density matrix p we have 
calculated the next quantities. Tracing out the matrix p over the field (atom) states we have 
calculated the reduced atom (field) density matrix p a ^ = Trf >a {p} , where Tr a j{p} are the 
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partial traces over the atom or field states accordingly. The latter allows us to study the 
dynamics of the excited-state probability P + (r) = p a ee , the mean photon number (n(r)) 
and the entropy S a (r) = —Tr{p a In p a }. Also the Fourier spectrum of P + {t) is explored. 
The convergence of the equations is tested and the dynamics of the system is studied for 
different values of the external field relative to the atom-field mode coupling. The results of 
the numerical solution of master equation (|19|) are shown in Figs^El 



III. NUMERICAL RESULTS 



We briefly consider the dynamics of a vacuum a = of driving JCM in order to un- 
derstand some general features first. FigC] shows the behavior of the atom quantities and 
mean photon number for resonant case (5 = 0) for vacuum initial state |0) at losses case, 
obtained as a result of the numerical solution of the master equation p9j) . This solution is 
close to Eq. (jl2)) . but with a damping due to both mirror losses in the cavity and sponta- 
neous emission from the atom. We use the following parameters £ =0.7, g = 0.2, 71/72 = 5 
and 71 = 5 • 10~ 3 . To make clearer the details of time evolution we use here a time inter- 
val r = 200 , which is large in comparison with the characteristic time scale of revival 25I 
r > 2t r , where r R = 2jn/g = 55, 7 = £/2g = 1.75. One can see, that dynamics of the 
probability of excited level occupations P + {j) has form of a damped sequence of collapses 
and revivals. Dynamical collapses and revivals are specific features of a unitary evolution. 
They are strongly affected by decoherency, which have the time constant set by the cavity 
field energy damping time ~ l/(7i + 72)- From Fig^a) one can see that both collapses and 
revivals are progressively less pronounced due to dissipation 7i )2 7^ 0. The dissipation also 
reduces the magnitudes of non-diagonal elements p a ge (in this case p a eg = (Pg e )* are purely 
imaginary quantities). One can see, that quantity Im(p^ e ) over time approaches to zero due 
to losses, which causes the field phase information to wash out and progresses decoherence 
in the coupled system. 

The dynamics Tr{(p a ) 2 } in Fig^a) shows that in the collapse area quantity Tr{(p a ) 2 } 
has a maximum close to 1 at r = t r /2 = 27.5. This point corresponds to the atomic 
attractor state which is completely independent on the initial atomic state [3]. At this 
point the compound system is in the disentangled state when in lossless system S a 'f — > 
and Tr{(p a 'f) 2 } — > 1. In lossy system such limit values are fulfilled only approximately. 
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One can see from FigOJb) that the mean photons number (n) increases from the initial 
zero value; however over time (n) assumes the steady-state value justified by the amplitude 
of external driving field £. Fig^c) shows the dynamics of the entropy of atom subsystem 
S a (r) = —Tr{p a In p a }. In the area of the first collapse S a has oscillating behavior. However 
these oscillations are smoothed away over time and 5' a (r) approaches to the steady-state 
value ln(2) = 0.69, which corresponds to maximally entangled state of two-level atom and 
field mode. At a lossless case the exact equalities Tr{(p a ) 2 } = Tr{(pf) 2 } and S a = in 
standard JCM hold [27[. Our simulations confirm this fact with very good accuracy. 

Next, we study details of the driving and initial coherent fields interaction. FigEf a) shows 
dynamics of the atom subsystem for the compensative case (7 = 7 + a = £/2g + a = 0) ( see 
Eq. (jT7Jl for the exact resonance case (5 = 0)), but taking into account dissipation 7^ 0. 
Comparison of Fig^a) and FigJ2(a) shows that the dynamics for this case is essentially 
different from the that in the initial vacuum case a = and has the form of damped vacuum 
Rabi oscillation in standard JCM. 

Dynamics of the mean photons number (n) is of interest and is shown in FigJ2Jb). The 
solid line in FigEfb) shows results (n) from a numerical simulation of the atom-field master 
equation (fT9"j) for the compensative case, taking into account losses in the system. For a short 
period of time r < 10 this simulation is in very good agreement with the exact formula f)18j) . 
However a discrepancy arises over time due to losses not included in f!18|) . One can see that 
despite of the mutual compensation of the initial coherent field and driving field in Eq. ([T8|) . 
the mean number of photons oscillates in vicinity ~ \a\ , that is justified by the field \S\. 
One can see from FigEfc), that for a short time the dynamics of entropy is similar to that 
in the vacuum case. However over longer time intervals the impact of dissipation becomes 
essential. Despite dissipation, at the time instances = n(2k + l)/2g , k = 0,1... the 
entropy S cl (ti : ) is very close to zero that implies the transition of the coupled atom-field 
system to the uncoupled pure state. At the moments Tk = nk/g the entropy has a value 
close to In 2, which corresponds to the maximum entanglement of the two-level atom and 
the field. Since the mean photon number does not vanish (see FigEfb)) this dynamics is 
stable asymptotically. 

FigEl shows dynamics of the coupled subsystem in the compensative case (7 = 0) when 
both dissipation 71,72 7^ and detuning 5^0 are non-zero. Notice in a non-resonant 
case (5 7^ 0) the time dependence in the Hamiltonian (J2J) is not removed, and therefore the 
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displacement operators transformation already does not allow to reduce the driving JCM to 
standard JCM. Compared to the case of the exact resonance (FigJ2J) the sequential collapses 
and revivals now practically disappear over long periods of time; in FigJ^a) the first collapse 
is seen only. Nevertheless some maxima of Tr{(p a ) 2 } remain resolvable. One can see from 
Figl^c) that during short periods of time r < 7 the behavior S a is similar to a that in the 
vacuum case (FigJ2fc)). However over longer periods of time the influence of decoherence 
becomes essential. The entropy quickly approaches its asymptotic value. 

Fig0] shows the Fourier spectra for cases of driving JCM considered above. This spec- 
trum, obtained by Fourier transform of numerically calculated P + (r), exhibits well separated 
discrete frequency components, which are scaled as square roots of the successive integers. 
The spectrum in FigEJ^a) corresponds to the time dynamics presented in Fig^a). One can 
see, that P + spectrum at the initial vacuum state |0) for driving JCM case is similar to that 
in a standard JCM (without driving field £ = 0) at initial coherent field \a). Note that for 
standard JCM such peaks were observed experimentally in Ref . j^fll] . The spectrum in Fig0] 
(a) is rather similar to the spectrum shown in Fig. 2(d) of experiment for 7 = 1.77. We 
can conclude that the vacuum case of driving JCM is similar to a standard JCM case with 
a coherent initial field, at least for the exact resonance. This conclusion reiterates the fact 
that coherent Held state is as dose to the classical fl e ld as q ua„tu m m echanics pe, m itsQ. 
Note that FigJ^a) shows that peak number 2 is the highest one. This also follows from the 
above mentioned. In this case 7 = 1.75 therefore the calculated number of highest peak is 
£[7 2 - 1] = 2. 

Figs|Ub)-|Hd) represent compensative case (fT7j) for the loss and detuning case. FigEfb) 
shows the spectrum for case 5 = 0, however for not small losses, which are. One can see 
only the peak, which corresponds to n = is present in the Rabi spectrum. Such spectrum 
in FigJH(b) corresponds to the measured spectrum shown in Fig. 2 (a) in Ref. |2o| for nearly 
vacuum case (no injected fields). Thus, the conclusion about the possibility subtracting of 
the coherent and classical fields remains valid in the loss case also. In FigEfc) influence of 
detuning is shown for smaller losses. One can see, that even for 5 = 0.1 which is not very 
small the spectrum of several first Rabi frequencies is well recognizable. There are a few 
peaks in FigE^c), which correspond to detuning 5^0. Such peaks are present even at a 
larger dissipation rate (the dissipation rate in Fig^d) are increased by two orders, with 
respect to Fig0Jc). Nevertheless the peak of vacuum oscillations remains dominant. 
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In general one can readily derive, that the quantity 7 can be rewritten (with the previously 

1 /2 

used notations) as 7 = S/2g = E e /E vac , where E vac = (frw a /2eoV) ' is the field per photon. 
Then from equality (fT7|) 7 = 7 + a = one can obtain E e = E va J\}l 2 exp(z(t> + 71")). The 
latter equation states that for the compensated case (jl7)l the driving field E e should be in a 
coherent state, but shifted by phase 7r concerning to initial state a. In general, condition (fTTjl 
does not hold exactly due to quantum fluctuations. However, since the ratio of fluctuation 
to the mean number (fractional uncertainty in the photon number, see e.g. Ch.3 in Ref. 3|) 
is An/n = n _1//2 for a Poisson process, the larger values of n become, the better is the 
accuracy of condition (fTTj) . 



IV. CONCLUSION 

Studied the field-atom interactions in the driving JCM shows that a variation of the 
amplitude or /and the phase of the driving field enables one to manipulate the dynamics 
and the spectrum of quantum Rabi oscillations. There are two distinct regimes. In case 
of summation of the driving field and initial coherent field one can underscoring a selected 
frequency in the Rabi frequency spectrum. The subtraction provides a possibility of com- 
pensation of both fields. For the case of the exact compensation, the frequency spectrum of 
a two-level atom becomes similar to that of vacuum oscillations in standard JCM. To use 
these processes for quantum information technology, the decoherence time must be greater 
than the time scale of the atom-field interaction. In this case, these processes may have vari- 
ous applications related to monitoring of the entanglement of two-level atom with quantized 
field and may be used to develop quantum information technology devices. 
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FIG. 1: Quantum dynamics in driving JCM as function of the interaction time r for initial vacuum 
case |0), S = A = 0, £ = 0.7, g = 0.2, 71 = 5 • 10" 3 , 72 = 10~ 3 . (a) Probability of the excited level 
|e) occupations P + (solid line), Im(/^ e ) (dashed line), Tr{(p a ) 2 } (dotted line), (b) Mean photon 
number (n). (c) Entropy S a (solid line), Tr{(p*) 2 } (dashed line). See details in text. 
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FIG. 2: Quantum dynamics in driving JCM as function of the interaction time r for initial coherent 
case | a), a = —1.75 (compensating case (|17)>). 5 = A = 0, £ = 0.7, g = 0.2, 71 = 5-10~ 3 , 72 = 10~ 3 . 
(a) Probability of the excited level |e) occupations P + (solid line), lm.(p^ e ) (dashed line), Re(/>^ e ) 
(dash-dot line), Tr{(p a ) 2 } (dotted line), (b) Mean photon number (n). (c) Entropy S a (solid line), 
Tr{(pf) 2 } (dashed line). See details in text. 
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FIG. 3: The same as in FigElbut for detuning 5 = 0.1. 
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FIG. 4: Fourier transforms of the probability P + (t) revealing the discrete Rabi frequencies, oc- 



curring at the successive square roots of the integers \/n + 1, numbers n = 0, 1,2... are inserted 
in vicinity of peaks. Spectra correspond to (a) Fig^a); (b) FigE(a); (c) FigEJa), but with 
7l = 5 • HT 5 , 72 = 10" 5 ; (d) FigUa). 
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